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In this paper, we investigate the dynamics of Born-Infeld(B-I) phantom model 
in the us — us' plane, which is defined by the equation of state parameter for 
the dark energy and its derivative with respect to 7V(the logarithm of the scale 
factor o). We find the scalar field equation of motion in us — u>' plane, and show 
mathematically the property of attractor solutions which correspond to us v ~ —1, 
Q. v = 1, which avoid the "Big rip" problem and meets the current observations 
well. 
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1. Introduction 



Recent observations of high-redshift Type la Supernova[l] and the Cosmic Microwave Background[2] have 
shown us such a fact: the density of clustered matter including cold dark matters plus baryons, Q m o ~ 1 /3, 
and that the Universe is flat to high precision, fltotai — 0.99 ± 0.03[3]. That is to say, we are living in a flat 
universe which it undergoing a phase of accelerated expansion, and there exists an unclumped form of energy 
density pervading the Universe. This unknown energy density which is called "dark energy" with negative 
pressure, contributes to two thirds of the total energy density. Perhaps the simplest explanation for these 
data is that the dark energy corresponds to a positive cosmological constant. However, the cosmological 
constant model suffers from two serious issues called " coincidence problem" and " fine-tuning problem" . An 
alternative is a scalar field which has not yet reached its ground state. These scalar field models include 
Quintessence[4-14], K-essence[15], Tachyon[16], Phantom[17-20], Quintom[21] and so on. 

On the other hand, the role of tachyon field field in string theory in cosmology has been widely studied[22]. 
It shows that the tachyon can be described by a B-I type Lagrangian resulting from sting theory. Nonlinear 
Born-Infeld scalar field theory, firstly proposed by W.Heisenberg in order to describe the process of meson 
multiple production connected with strong field regime [23], is discussed in cosmology recently [24]. It shows 
that the lagrangian density of this B-I scalar field posses some interesting characteristics[25,26]. In Ref[25], 
the author showed that a horizon exists for a large class of solution in which the scalar field is finite though this 
horizon is singular. Naked singularities with everywhere well-behaved scalar field in another class of solution 
have also been found. Lately the quantum cosmology with the B-I scalar field has been considered [27], in 
the extreme limits of small and large cosmological scale factor the wave function of the universe was found 
by applying the methods developed by Vilenkin, Hartle and Hawking. The result has suggested a non-zero 
positive cosmological constant with largest probability, which is consistent with current observational data. 
The classical wormhole solution and wormhole wavefunction with the B-I scalar filed has been obtained in 
Ref[28]. The phantom cosmology based on B-I scalar field with a special potential has been considered in 
Rcf[29]. The result shows that the universe will evolve to a de-sitter like attractor regime in the future 
and the phantom B-I scalar field can survive till today without interfering with the nucleosynthesis of the 
standard model. 

Recently, many authors have investigated the evolution of dark energy models in the us — u>' plane[30,31], 
where us' is the time variation of u> with respect to N . According to different regions in the us — us' phase 
plane, these models can be classified two types which are call "thawing" and "freezing" models. Based 
on this consideration, we study the cosmological dynamics of B-I type phantom model with the potential 
u(ip) = Aipe~^ v in the ui — us' plane and investigate the properties of attractor solutions. According to the 
obtained B-I field equation of motion, we analyze the critical points and find the critical point with u> ~ — 1 
is the late-time attractor, where B-I field becomes ultimately frozen, as shown in the Fig.l. 



2. B-I Type Phantom Model 

In 1934[32], Born and Infeld put forward a theory of non-linear electromagnetic field. The lagrangian density 
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The lagrangian density for a B-I type scalar field is 
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Eq.(2) is equivalent to the tochyon lagrangian [— V{ip)sjl — g^f. ^(p, v + A] if V(ip) = — and cosmological 

constant A = — (— is two times as "critical" kinetic energy of if field). The lagrangian (2) possesses some 

interesting characteristics, it is exceptional in the sense that shock waves do not develop under smooth 
or continuous initial conditions and because nonsingular scalar field solution can be generated [33]. When 
rj — > 0, by Taylor expansion, Eq.(2) approximates to the lagrangian of linear scalar field. 
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Now we consider the Lagrangian of phantom model with a potential u(ip) in spatially homogeneous scalar 
field, Eq.(2) becomes 
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In the spatially flat Robertson- Walker metric ds 2 = dt 2 — a 2 (t)(dx 2 + d 2 y + d 2 z), Einstein equation = 
KT^ V in phantom dominated epoch, can be written as 
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where p v and p v are the effective energy density and effective pressure respectively, and we work in units 
8ttG = 1. The Energy- moment tensor is 
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From Eq.(7), we have 
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Pv = lu vPv = —T{ = -T 2 = -T* = 
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From Eqs.(8)(9), we get 
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It is clear that the equation of static lo < — 1 is completely confirmed by Eq.(10)and it accords with the 
recent analysis of observation data. 



3. Dynamics of B-I Type Phantom Model In w - w' Plane 

We make such a definition TV = Ina which results in the no-dimension of scale factor a. From Eqs.(5)(6)(9) 
and (10), we have 
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Now, we defining a function 
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So, we can rewrite Eq.(14) as 
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and the sign ""' denotes the derivative of u(ip) with respect to if. Eq.(15) is the scalar field equation of 
motion in B-I type phantom model. For the B-I field rolling down its potential, the ± sign before u'((p) 
corresponds to u 1 {(f) > or u'((p) < 0. When 77 — > 0, the scalar field equation of motion in B-I type phantom 
model will tend to be the canonical case. 

Now we consider the potential u(ip) = A(pe~^ v . The system decided by Eq.(15) admits these critical 
points lu v = 1, lu v ~ — 1, A = 1. When A = 1, uj a varies very slowly where the B-I field is tracking and the 
ratio of kinetic energy to potential energy of the B-I field becomes a constant. When uj^ ~ —1, it corresponds 
to a late time attractor where the B-I field becomes ultimately frozen, as shown in the Fig.l. When we set 
different initial conditions, the critical point always tends Lo a ~ —1, as shown in Fig. 2. Fig. 3 shows the 
evolution of the state equation of state uj v oscillates at the beginning and tends to be -1 with respect to N 
eventually. 




Fig.l The attractive behavior of B-I model with potential 
u((p) = A(pe~ ,3 ' p in the u — uj' plane. We set rj — 10~ 4 , 
13 = 0.5 and A = 1.0. 
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Fig. 2 The comparison of attractive behavior of B-I 
model with potential u(ip) = ALpe~^ v for different initial 
conditions in the u> — u' plane. 
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Fig. 3 The evolutive behavior of equation of state parame- 
ter Uip. Uip — 1 which corresponds to a late time attractor, 
where the B-I phantom field becomes ultimately frozen. 
| 1 + |=| 1+ (-1.08) |= 0.08 > 0.02, which satisfies 
the limit of "freezing" model. 



4. Conclusions 

Caldwell and Linder[29] classified the quintessence models into two types "thawing" and "freezing" model 
according to the different regions in ui — lo' plane and gave the limit of quintessence [28]: 1 + u> > 0.004 for 
"thawing" model and 1 + lo > 0.01. We consider the dynamics of B-I phantom model in the lo — lo' plane for 
the first time. In this paper, we investigate the cosmological dynamics of B-I phantom model with potential 
ALpe~^ v in the lo — lo' plane and examine the evolution of equation of state parameter with respect to N. 
These numerical results show that the critical point with lo ~ — 1 is a late-time attractor, where B-I field 
becomes ultimately frozen and lo v ~ — 1. The existence of the attractor solution prevents the phantom 
energy from increasing up to infinite in a finite time, therefor the presence of phantom energy dose not lead 
to a cosmic doomsday in a theory with late time attractor. We also consider the effect of the initial condition 
on the property of attractor solution. | 1 + lo v | = | 1 + (—1.08) |= 0.08 > 0.02, which satisfies the limit of 
"freezing" model. So, we know our B-I model belongs to "freezing" model, where the B-I field which was 
already rolling up towards its potential maximum -§g which is decided by u'(tp) — 0, prior to the onset of 
acceleration, but which slows up and creeps to a halt as it comes to dominate the Universe. 
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